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In this paper, we studied the doped Z2 topological spin liquid of the toric-code model. We found 
that the doped holes become supersymmetric particles. The ground state of the doped Z2 topological 
spin liquid becomes new matters of quantum states - supersymmetric Bose-Einstein condensation 
or supersymmetric superfluid. As a result, this system provides a unique example of manipulatable 
supersymmetric many-body system. 



Doped holes in a spin model had become an impor- 
tant issue since the discovery of high Tc superconduc- 
tivity in cupratesQ. As the microscopic model of the 
high Tc superconductivity - the t — J model has been 
intensively studied for several decades. Motivated by 
the experimental facts in the high-Tc cuprates, the spin- 
charge separation idea [3] was a very basic concept by 
introducing spinless "holon" of charge e and neutral spin- 
1/2 "spinon" as the essential building blocks of the re- 
stricted Hilbcrt space. However, for a long range anti- 
ferromagnetic (AF) order, there is no true spin-charge 
separation at all. People pointed out that a single doped 
hole in the AF order can be a charged spin bagQ, or 
Shraiman-Siggia dipolejl], or localized objectfl, Q from 
different points of view. The true spin-charge separa- 
tion occurs only in the quantum disordered spin states 
(people call them quantum spin liquid states). And the 
doped spin liquid is always a superconducting order with 
holon-condensation. 

On the other hand, in the last decade, several ex- 
actly solvable spin models with Z2 topological spin liquid 
were found, such as the toric-code model 0, the Wen- 
plaquette model [1, Q and the Kitaev model on a hon- 
eycomb lattice [l^l- It becomes an interesting issue to 
study the properties of doped spin liquid by doping holes 
to these exactly solvable spin models. In Ref . . the 
fermi liquid nature of doped Z2 topological spin liquid is 
obtained. People also studied the effect of doped holes in 
the gapless phase of the Kitaev model in Ref. and 
pointed out that the topological superconducting state 
can be its ground state. 

In this paper, after studying the non-perturbative 
properties of doped holes, we found that a doped hole 
in the Z2 topological spin liquid of the toric-code model 
becomes a supersymmetric holon and a universal feature 
of the doped Z2 topological spin liquid is the emergence 
of super symmetry. Wc found that the ground state of the 
doped Z2 topological spin liquid is supersymmetric Bose- 
Einstein condensation (SBEC) state, of which there exist 
the fermionic Goldstone mode - Goldstino which is the 
partner of the Bosonic Goldstone mode. And we can ma- 



nipulate the supcrsymmetry by tuning the transverse ex- 
ternal field to the system. After breaking the supcrsym- 
metry by the transverse external field, the ground state 
becomes a new matter of a quantum state - supersymmet- 
ric superfluid - Bose-Einstein condensation (BEG) state 
for the bosonic holons and superconducting (SC) state for 
the fermionic holons. That is, the doped Z2 topological 
spin liquid provides a unique example of supersymmetric 
many-body system. 

Our starting point is the so-called t-toric-code model 
that is described by the following Hamiltonian 



H — Ht + Htc + Hi , 



(1) 
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where Zi 

X,=«+e.<+e.+e„<+e„ with A > B > Q. a 

are the spin-indices representing spin-up (cr =t) and 
spin-down (ct =4,) for the electrons. ^ is the chemical 
potential, (i, j) denote two sites on the nearest-neighbor 
link. Vs projects the Hilbert space onto the subspace 
of one electron per-site. s^'*''^ are the spin operators 
of the electrons. For simplify, we consider the case of 
A = B = g in this paper. Hi is the external field terms 
with /i^ > 0, hy > 0. 

The ground state of Htc denoted by Zi = +1 and 
= +1 at each site is a Z2 topological state. For 
the toric-code model, the elementary excitations are Z2 
vortex {Zi = — 1 at even sub-plaquettc) and Z2 charge 
(Xi = — 1 at odd sub-plaquctte) . Fcrmions arc the bound 
states of a pair of Z2 vortex and Z2 charge. When we 
add the transverse external field, the total spin Hamil- 
tonian Htc + Hi cannot be solved exactly. However, for 
small value of the external field, the ground state is still 
Z2 topological order and the topological properties don't 
change. 

Internal degree of freedom of a hole and emergent su- 
pcrsymmetry : Now we study a single hole doped in the 
toric-code model. As shown in Fig.l, a hole is a charged 



2 





Hollow region 



FIG. 1: The illustration of a single hole. 



vacancy to remove the spin degree of freedom at site i, as 
c^j — >■ c|. Then near a hole, four plaquette operators are 
forced to be zero, i.e. .^i-e, = ^j-e„ = Xi = X^.g^-e^ = 
0. Around this hole, there exist a boundary that sepa- 
rates the topological spin liquid outside and a hollow area 
around the hole. For a hole in the Z2 topological order, 
the ground state degeneracy becomes 2 (without being 
considered the edge states) [l^ [i3 • To classify the degen- 
eracy of the ground states, as shown in Fig. 2, we define 
two types of closed string operators, Wy{CA) = Oc *f 
and W}{Cb) = Oc^f- ^v{'^a) is a closed string opera- 
tor around a hole. W/(C_b) is a closed string operator 
from the boundary of a hole to the boundary of the sys- 
tem. Since Wi,(Ca) and Wf{CB) form the Hciscnberg 
algebra, we can map the Hilbert space of Wv{Ca) and 
H^/(C_b) onto that of a pseudo-spin [S = \) and iden- 
tify WyiCA) and VF/(Cb) to spin- 1/2 operator and 
as WJCa) 5^ and WfiCs) ^ 5^. Then, we use |t) 
and 14.) to denote the pseudo-spin degree of freedom of 
the quantum states of a hole. These two quantum states 
can be characterized by the fermion parity (or the flux 
quanta of the Z2 vortex) of the hole : |t) corresponds to 
the state with even fermion parity, \\.) corresponds to the 
state with odd fermion parity. Now the hole's statistics 
dependents on its fermion parity. 

For the zero external field case, — — 0, the two 
quantum states |t) and H) are degenerate exactly. When 
wc add the external field, the quantum tunneling effect 
will lead to an energy splitting. Then wc can derive the 
effective model of a hole by the quantum tunneling theory 

as Tih — h^S^ + h^S^ where = ^^gg"L°li and = 

^-4g)''ii ■ Here is the distance between the boundary 
of the hole and the boundary of the system. See the 
detailed calculations in Ref.fli, [l^. In general, for an 
infinite system, we have = 0. In addition, for the 
two-hole case, we have the effective model of two holes 




string W,(C_) hole 

FIG. 2: The illustration of holes in a Z2 topological order of 
the toric-code model. 



as 

=^ .F'^S'^S^ + r-'^SlS'i + hfS^ (2) 

1=1,2 

where J^^ = ^I'sjy-i , J^'' - jz^ji^TT^- I is the distance 
between the holes and L^^ is the distance for a path sur- 
rounding the two holes. Due to J^^ ^ J^^, the effective 
model of the two-hole case is reduced into 

-Hir ^ J'^'-'SfS'^ + h'-Sl + h'^SI . (3) 

Due to the term of J^^5fiSf , two holes may exchange 
their fermion parities. 

So one may tune each parameter in above effective 
Hamiltonian Hh by controlling the external field along 
special direction and then manipulate the internal degree 
of freedom of holes. For example, in the limit J^^ — )■ 0, 
the quantum state with lowest energy is |4,)2 
with odd fermion parity for each hole. In the limit, 
0, the eigenstates are (|t)i+ |4-)i) (lt)2+ |l)2) 
and (It)i- |i)i)®(|t)2- |4-)2)- Now the average fermion 
parity of each hole is 1/2 for both eigenstates. 

Effective super symmetric model: Because the original 
statistics of a hole is fermion, a holon (a hole in the Z2 
topological state) is a boson for the internal state \\^) with 
odd fermion parity or a fermion for the internal state 
It) with even fermion parity. Then for a single holon, 
we can use a two-component supersymmetric operator to 
describe it as c^i. c^i plays the role of a two-component 
spinor with a pseudo-spin where 0:^.^ = bi and cj,- = fi are 
the charged boson and charged fermion operators on site 
i. 

Now we consider the t-toric-code model with finite 
hole-concentration, of which the effective Hamiltonian is 
Hcs = Ht + Hh where Ht is the hopping term of super- 
symmetric holons, 

Ht = E 'Pd^^-^J^'s+^iB E blh+flF J2 flf^- (4) 
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Vs projects the Hilbert space onto the subspace of one 
particle per site. Now the projector Vs has no effect 
for the spinless fcrmionic holon while guarantees the 
particle number of bosonic holon n\ = b\bi. /jb and fxp 
are the chemical potentials for the bosonic holons and 
that of fermionic holon, respectively. 

To characterize the supersymmetry, we introduce the 
generators Sf , Sf , Si that are just the string opera- 
tors WfiCs), -iWf{CB)Wv{CA), Wy{CA)- And we have 
5f = {S-+S+)/2, Sf = {S--S+)/2i. It is clear that S^ 

is a fermionic operator. Physically, iS~ turns a fermionic 
holon into a bosonic holon, and S^ does the opposite. By 
the two-component supersymmetric operator, we have 
S- = h\U. St = hfl and St = (blh - flh)/2. Thus 
we re-write the term T^h into 

ij i 

where J^^ = jz^^y^- Because the exchange term 
decays exponentially, we may only consider the shortest 
case that is / = 2 (for I < 2, the hollow regions of two 
holons merge and we cannot define the corresponding 
quantum tunneling process). As a result, the fermion 
parity exchange term J^ij Ji '^SfS^ is reduced into 

5r+3,^ + 5f 5f+3e„ ) (6) 

i 

i 

+ {bl.h + hfI){blse,.U+3e, + 

Thus we have an exact supersymmetry for the toric 
code model without the external field, i.e., [S~ , i?cff] = 
[iSj^, Hcff] = 0. When there exists the external field {h^ ^ 
0, ^ 0), the supersymmetry is broken explicitly, i.e., 
[Sr, H,ff]^0, i/eff] ^0. 

Supersymmetric Bose-Einstein condensation: We 
firstly study the many-body system with exact super- 
symmetry as [S^, H] = 0. Now for the dilute gas limit 
of holons (the hole-concentration 5 is smaller than 1%), 
the effective model is 

HeS = -t E bib, - t ^ ///, +^,bY1 blb^ + l^F flh 
{ij) (JJ> » » 

(7) 

where fip = fJ'B = — 4t. bj denotes the annihilation oper- 
ator of the hard-core bosons. Because in the small hole- 
concentration limit, we can release the single-occupied 
condition of the bosonic holons and consider the holons 
as non-interacting bosons. For non-interacting bosons 
we always have fiB = — 4t (we measure the energy from 
the bottom of single particle dispersion) at zero tem- 
perature, regardless of boson number. Supersymmetry 
of iJoff requires fip = fiB, as a result we have or 1 
fcrmionic holon. Now each holon becomes a boson and 



the ground state is a BEC state with holon-condensation 
as (0|fek=o|0) = &oe"^° ^ 0. Here |0) denotes the ground 
state. Thus the ground state spontaneously breaks both 
global U{1) symmetry and supersymmetry. 

There arc two types of collective excitations. One type 
is the Goldstone mode that describes the phase fluctua- 
tions as ipQ ^ Lpo + fi- The effective Hamiltonian of the 
Goldstone mode is H = — 2<&oS(y) cos(v3i — ipj) which 
is really a two-dimensional XY model. Another type of 
collective excitation is the Goldstino mode which can be 
regarded as the fermionic " spin wave" [l^l . Due to the 
holon-condensation, we have S^ = bofj. From the hop- 
ping term of // , we have the effective model of the Gold- 
stino as i?Goidstino = -t\bo\~'^J2{ij)'Si''Si- The disper- 
sion of the Goldstino mode is given by ~ \bo\ ^ e-^ 
where £k = — 2i(cos + cos ky). Due to the existence of 
the Goldstino mode, we call the unique ground state to 
be supersymmetric Bose-Einstein condensation (SBEC). 

Supersymmetric Bose-Einstein condensation with 
Fermi surface: For the case of 0, = 0, the 

effective model turns into 

^^eff = -t E - * E ///j- + E ^Bbjb, 

(ij) (ij) i 

+ Ea'^^///^ + yE(^'^»-///'0 (8) 

i i 

where = -^^Tj^yrr- The supersymmetry is partially 

broken by term. 

Because there is no couple between fermionic holons 
and bosonic holons, the ground state can be considered to 
be a Bosc-Fcrmi mixture. The bosonic holons condense 
and fermionic holons form Fermi liquid. In the dilute hole 
limit, we can simplify the Fermi surface of the fermionic 
holons to be a circle with a radius kp- Now in continuum 
limit, the effective Hamiltonian can be reduced into 

^ E[£ + if^BUbib^ + E[£ + i^^^Uflh 

k k 

(9) 

where the effective chemical potential of bosonic holons 
{fiB)cs and the effective chemical potential of fermionic 

holons (/^F)eff are i^J.B)es = Ms + ^ and {^J.F)cS = - 
respectively, m is the effective mass of holons, m ~ 

j_ 

2t ■ 

Then we can estimate the number of fermionic holons 
and that of bosonic holons by minimizing the total 
ground state energy. We define the total holon number 
to be Nt = Nb + Nf where Nb = J2z and Nj = ■ 
The total ground state energy is i?totai = E^ + Ef, where 
Eb ~ ^Nb ~ ^{Nt — Nf) is the ground state energy of 

bosonic holons and Ef = ((5/)^ ^ — — y~ is the ground 
state energy of fermionic holons. Here N is the lattice 
number and ^/ is the concentration of fermionic holons. 
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completely. Now the effective model is 
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FIG. 3: The phase diagram: the supersymmetric Bose- 
Einstein-condensation (SBEC) (the blue line at = 0), the 
supersymmetric Bose-Einstein-condensation with Fermi sur- 
face (SBEC+FL), the Fermi liquid (FL). There is a quantum 
phase transition at {h^)c (the red line). 



H = -tY^ bib, - t J2 flh + A'B E 



+ 



J; 



+ {b\h + bd}){bl,.f 



(10) 



The ground state is really a BEC of bosonic holons 
(0|5o|0) = boe'-'^'' with induced pairing of fermionic 
holons. 

Let us estimate the induced SC pairing of fermionic 
holons. The effective Hamiltonian of fermionic holons is 



(11) 



From the condition, ^gj^*'^' = 0, we obtain the concentra- 
tion of fermionic holons as Sf = ^^^^ = Because the 
maximum fermionic holon's concentration is i5 which is 
the hole concentration, we get a critical value of the ex- 
ternal field, {h')c = 'iTTtd or (/i^)^ = (4<?)"/^^ {^y^^^- 

For small external field case, < {h^)c, there exists 
bosonic holons with the concentration to he Sb = S — S j . 
Now the ground state is a mixture of SBEC of bosonic 
holons and Fermi liquid (FL) of fermionic holons (we call 
it SBEC-I-FL state). In general, is very tiny. For 
example, for the case of /g = 0.3, we have /g — S x 
10~^^. So we always have the SBEC state for bosonic 
holons and the Fermi liquid for fermionic holons with 
very tiny fermi surface. 

For large external field case, > {h'^)c, all holons 
are fermions. The ground state is a Fermi liquid (FL) 
state with the number of fermionic holons as TV/ = Nt- 
Now the effective chemical potential of fermionic holons 
is {^p)^ff ~ 2TTt6 and the effective chemical potential of 
bosonic holons is {fiB)cff — ^^/2. There exists a finite 
energy gap to excite a bosonic holon. While the fermions 
have no energy gap. 

From these results we plot the phase diagram in Fig. 3, 
of which the red line denotes the quantum phase tran- 
sition between SBEC-I-FL and FL. The FL state is only 
stable in the limit of t(5 — and /g — > 1. However, 
for > 0.34(7, a topological quantum phase transition 
occurs and the ground state turns into a spin polarized 
state without topological order [T^423j . Thus we can only 
discuss the case of small external field as < 0.34(7. 

Supersymmetric superfluid: In this section we study 
the case of 0, 0. The supersymmetry is broken 



Thus, when the bosonic holons condensate, there exists 
Px+ Py superconducting order parameter for the non- 
interacting fermionic holon 



(0|/,I/lk|0> 



l2 txx 
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[(sin 3^2; + sin 3ky)] . 



(12) 



Now we have a gapless fermionic holon with Pa;-|- Pj, pair- 
ing. 

Due to the condensation of the bosonic holons and 
the induced SC pairing of fermionic holons, the effective 
model of the Goldstino mode turns into 



Goldstino 



(13) 



Then the dispersion of the Goldstino mode is derive as 
Ek = l^oP" v^e| + A| where Sk — ~2t{cos k^^ + cos ky) + 
{^lF)^s and A2 = (ffl^)2[(sin3fc:, -hsin3fcy)2]. 

Conclusion : In this paper, we found that for the 
Z2 topological spin liquid of the toric-code model, a 
doped hole becomes a charged supersymmetric parti- 
cle. Thus the doped Z2 topological spin liquid of the 
toric-code model provides a unique example of super- 
symmetric many-body system. The ground state is a 
new matter of quantum state - SBEC, of which there ex- 
ists fermionic Goldstone mode - Goldstino which is the 
partner of Bosonic Goldstone mode. And we can tune 
the supersymmetry by adding transverse external field 
to the toric-code model. After breaking the supersym- 
metry by transverse external field, the ground state may 
be a supersymmetric superfluid - BEC state for bosonic 
holons and SC state for fermionic holons. 
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